Dynamical creation of a supersolid in asymmetric mixtures of bosons 
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We propose a scheme to dynamically create a supersolid state in an optical lattice, using an attrac- 
tive mixture of mass-imbalanced bosons. Starting from a "molecular" quantum crystal, supersolidity 
is induced dynamically as an out-of-equilibrium state. When neighboring molecular wavef unctions 
overlap, both bosonic species simultaneously exhibit quasi-condensation and long-range solid order, 
which is stabilized by their mass imbalance. Supersolidity appears in a perfect one-dimensional 
crystal, without the requirement of doping. Our model can be realized in present experiments with 
bosonic mixtures that feature simple on-site interactions, clearing the path to the observation of 
supersolidity. 
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The intriguing possibility of creating a quantum hy- 
brid exhibiting both superflow and solidity has been en- 
visioned long ago [1 . However, its experimental observa- 
tion remains elusive. The quest for supersolidity has been 
strongly revitalized by recent experiments showing possi- 
ble evidence for a non-zero superfluid fraction present in 
solid ^He [2 . Yet, several theoretical results [3] appear 
to rule out the presence of condensation in the pure solid 
phase of ^He, and various experiments [4 show indeed a 
strong dependence of the superfluid fraction on extrinsic 
effects, such as ^He impurities and dislocations. While 
the experimental findings on bulk ^He remain controver- 
sial, optical lattice setups [5 offer the advantages of high 
sample purity and experimental control to directly pin 
down a supersolid state via standard measurement tech- 
niques. A variety of lattice boson models with strong 
finite-range interactions has been recently shown to dis- 
play crystalline order and supersolidity upon doping the 
crystal state away from commensurate filling |3 ; yet 
sizable interactions with a finite range are generally not 
available in current cold-atom experiments. Such inter- 
actions can be in principle obtained effectively by adding 
a second atomic species of fermions [7 , which, however, 
does not participate in the condensate state, in a way 
similar to the nuclei forming the lattice of a supercon- 
ductor without participating in the condensate of elec- 
tron pairs. 

Here we demonstrate theoretically a new route to su- 
persolidity, realized as the out-of-equilibrium state of 
a realistic lattice-boson model after a so-called "quan- 
tum quench" (a sudden change in the Hamiltonian) . 
The equilibrium Hamiltonian of the model before the 
quench realizes a "molecular crystal" phase character- 
ized by the crystallization of atomic trimers made of two 
mass-imbalanced bosonic species. Starting from a solid 
of tightly-bound trimers and suddenly changing the sys- 
tem Hamiltonian, the evolution induces broadening and 
overlap of neighboring molecular wavefunctions leading 
to quasi-condensation of all atomic species, while crys- 
talline order is maintained (Fig. [T]). Our model requires 
only local on-site interactions as currently featured by 
neutral cold atoms, which make the observation of a su- 
persolid state a realistic and viable goal. 

We consider two bosonic species {a =T, i) tightly con- 
fined in two transverse spatial dimensions and loaded in 
an optical lattice potential in the third dimension. In the 
limit of a deep optical lattice, the dynamics of the atoms 
can be described by a model of lattice hardcore bosons 
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FIG. 1: Dynamical onset of supersolidity by quantum 
quenching a mixture of light and heavy bosons, (a) A 

product state of bosonic trimers is the initial state of the evo- 
lution (larger symbols represent the |-bosons); switching off 
one of the superlattice components leads to a supersolid state 
in which the particles delocalize into a (quasi-) condensate 
while maintaining the original solid pattern without imperfec- 
tions, (b) Momentum profile of the |-bosons, (n|.) vs. time 
in units of hopping events h/ J^. A quasi-condensate peak de- 
velops rapidly. Inset: Density distribution (nj) averaged over 
the last third of the evolution time, showing that crystalline 
order is conserved in the system. (The simulation parameters 
are L = 28, = 18, A^t = 9, J|/Jt = 0.1, U/J^ = 3.0.) 

in one dimension ^ |9] 

H = - ^ J. + h.c.) -UY^ n,,rn,,|. (1) 

i,cr i 

Here the operator b]^ {hicr) creates (annihilates) a hard- 
core boson of species a on site z of a chain of length 
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FIG. 2: Phase diagrams in and out of equilibrium. 

(a) Equilibrium phase diagram (empty circles). The dash- 
dotted line represents the points where the hopping of the 
j-bosons, J|, overcomes the energy gap to crystal disloca- 
tions, giving rise to the solid/super- Tonks (s-Tonks) transi- 
tion. The dashed line marks the points where a single-trimer 
wavefunction spreads over 2.8 sites, (b) Out-of-equilibrium 
phase diagram. An extended supersolid phase exists in the 
transient state attained after the quantum quench. Blue sym- 
bols delimit the boundaries of the solid phase, red symbols 
mark the lower boundary for the quasi-condensed (q-c) phase. 
The overlap of both phases (blue shaded region) is identi- 
fied as the supersolid phase. The yellow- filled symbols cor- 
respond to equilibrium data points. The lower boundary of 
the superfluid/super- Tonks region of the equilibrium phase 
diagram is seen to coincide with the lower boundary of the 
supersolid region out of equilibrium. 

L, and it obeys the on-site anticommutation relations 

{6icr,^Jcr} — 1- ^*cr = ^\(j^ia IS the numbei operator. 
Throughout this Letter we restrict ourselves to the case 
of attractive on-site interactions /7 > and to the case 
of mass imbalance, J| > J|. Moreover we fix the lattice 
fillings of the two species to = 1/3 and = 2/3. 

In the extreme limit of mass imbalance, J| = 0, ([T]) re- 
duces to the well-known Falicov-Kimball model of mobile 
particles in a potential created by static impurities [10^. 
For the considered filling it can be shown via exact diago- 
nalization that, at sufficiently low attraction U / < 2.3, 
the ground state realizes a crystal of trimers formed by 
two |-bosons "glued" together by an |-boson in an atomic 
analogue of a covalent bond (see Fig. [l^ for a scheme 
of the spatial arrangement). The trimer crystal is pro- 
tected by a finite energy gap against dislocations of the 
|-bosons, and hence it is expected to survive the presence 
of a small hopping J|. We have tested this hypothesis 
with extensive quantum Monte Carlo simulations based 
on the canonical Stochastic Series Expansion algorithm 
[TTJ [12]. Simulations have been performed on chains of 



size L = 30, ...,120 with periodic boundary conditions, 
at an inverse temperature pj^ = 2L/3 ensuring that the 
obtained data describe the zero-temperature behaviour 
for both atomic species. 

Fig. |2^ shows the resulting ground-state phase dia- 
gram, which features indeed an extended trimer crystal 
phase. For U/J^ > 2.3, and over a large region of JJJ^ 
ratios, the ground state shows instead the progressive 
merger of the trimers into hexamers, dodecamers, and 
finally into a fully collapsed phase with phase separation 
of the system into particle-rich and particle- free regions. 

For U/J^ ^ 2.3, increasing the Ji/J] ratio allows to 
continuously tune the zero-point quantum fluctuations of 
the |-atoms in the trimer crystal and to increase the ef- 
fective size of the trimers, whose wavefunctions start to 
overlap. We find that, when trimers spread over a criti- 
cal size of 2.8 lattice sites, they start exchanging atoms 
and the quantum melting of the crystal is realized. The 
melting point is also consistent with the point at which 
the hopping J| overcomes the energy gap to dislocations 
(dash-dotted line in Fig. |2|l). The resulting phase af- 
ter quantum melting is a one-dimensional superfluid for 
both atomic species: in this phase quasi-condensation 
appears, in the form of power-law decaying phase corre- 
lations {h\ ^hj^cr) oc \ri — Tjl""'^, which is the strongest 
form of oflF-diagonal correlations possible in interacting 
one-dimensional quantum models [13^ . Yet in the super- 
fluid phase strong power-law density correlations survive, 
(^i,cr^j,cr) 0^ cos(g'tr(n — ^j)) \Ti—rj\~^''^ exhibiting oscil- 
lations at the trimer-crystal wavevector q^^ = 27r/3. Such 
correlations stand as remnants of the solid phase, and in 
a narrow parameter region they even lead to a divergent 
peak in the density structure factor, S'cr(^tr) c<: L^"" with 
< /^cr < 1, where 

This phase, termed "super- Tonks" phase in the literature 
on one-dimensional quantum systems [14 , is a form of 
quasi-supersolid, in which one-dimensional superfluidity 
coexists with quasi-solid order. (Notice that true solidity 
corresponds to /3cr — 1.) 

The strong competition between solid order and su- 
perfluidity in the ground-state properties of this model 
suggests the intriguing possibility that true supersolidity 
might appear by perturbing the system out of the above 
equilibrium state. In particular we investigate the Hamil- 
tonian evolution of the system after its state is prepared 
out of equilibrium in a perfect trimer crystal. The initial 
state is a simple factorized state of perfect trimers (see 
Fig.[l}i): 

L/3 

l*o) = (g)|*i?"-^^) (3) 

n=l 

where the trimer wavefunction reads 

This state can be realized with the current technology 
of optical superlattices [15 , by applying a strong second 
standing wave component V^^ cos^[(/c/3)x + 7r/2] to the 
primary wave, V^-^ cos^(/cx), creating the optical lattice 
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along the x direction of the chains. This superlattice po- 
tential has the structure of a succession of double wells 
separated by an intermediate, high-energy site. Hence 
tunneling out of the double wells is strongly suppressed, 
stabilizing the factorized state (|3|. After preparation 
of the system in the initial state, the second compo- 
nent of the superlattice potential is suddenly switched 
off {Vx2 0) and the state is let to evolve with the 
Hamiltonian corresponding to different parameter sets 
([// J| , J| / J| ) . The successive time evolution over a short 
time interval [0, r] with r = 3h/ is computed using 
the Matrix-Product-States (MPS) algorithm on a one- 
dimensional lattice with up to 28 sites and open bound- 
ary conditions [16]. A bond dimension D = 500 ensures 
that the weight of the discarded Hilbert space is < 10~^. 
The evolution time step dt = 5x 10~^h/ is chosen so as 
to make the Trotter error smaller than 10~^. We charac- 
terize the evolved state by averaging the most significant 
observables over the last portion of the time evolution 
r/3. 
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FIG. 3: Coexistence of solid order and quasi- 
condensation in the supersolid phase, (a) The struc- 
ture factor peak *S'(^tr = 27r/3) scales linearly with system 
size L, demonstrating solid order for both bosonic species, 
(b) The density peak in momentum space (tiI^q) is plot- 
ted vs. L on a log-log scale, showing algebraic scaling and 
thus quasi-condensation. Boxes (diamonds) stand for par- 
ticle species [ (t), respectively. The data represented by 
blue boxes in part (a) is offset by -0.2 for better visibility. 
Parameters: J|/Jt — 0.1, [//J| = 3.0 (blue symbols) and 
J|/J| = 0.15,?7/J| = 2.5 (red symbols), (c) Square mod- 
ulus of the natural orbital Xi^"^ corresponding to the largest 
eigenvalue of the OBDM, calculated at final time r. In the 
supersolid regime (blue/red symbols for |/t bosons), the nat- 
ural orbital shows the characteristic crystalline order. This 
pattern is washed out in the purely quasi-condensed regime 
(dashed/solid curves for i/T)- The supersolid data is offset 
by +0.02 for the sake of visibility. Parameters: Ji/J^ — 0.1 
(supersolid), J|/Jt — 0.8 (quasi-condensed), U/J^ = 3.0, 
= 18, =9, L = 28. 

We find three fundamentally different evolved 
states, whose extent in parameter space is shown 
on the non-equilibrium phase diagram of Fig. ^p: 
Firstly, a superfluid phase, in which the initial crys- 
tal structure is completely melted by the Hamiltonian 
evolution, and coherence builds up in the system lead- 
ing to quasi-condensation out-of-equilibrium, namely to 



the appearence of a (sub-linearly) diverging peak in the 
momentum distribution (n^) = e^^^'^^~^^\h\^^hj^a) 
at zero quasimomentum, (ti^^q) oc L"'^ with < 
aa- < 1. Despite the short time evolution, quasi- 
condensation of the slow |-bosons is probably assisted 
by their interaction with the faster |-bosons, and 
is observed to occur for all system sizes considered. 
Secondly, we find a solid phase, in which the long-range 
crystalline phase of the initial state is preserved, as 
shown by the structure factor which has a linearly di- 
verging peak at the trimer-crystal wavevector S{qtr) oc L. 
Thirdly, an extended supersolid phase emerges, with per- 
fect coexistence of the two above forms of order for both 
atomic species. This is demonstrated in Fig. |3] a,b via 
the finite-size scaling of the peaks in the momentum dis- 
tribution and in the density structure factor. In this 
phase, which has no equilibrium counterpart, the Hamil- 
tonian evolution leads to the delocalization of a signifi- 
cant fraction of j- and |-bosons over the entire system 
size. Consequently quasi-long-range coherence builds up 
and the momentum distribution, which is completely flat 
in the initial localized trimer-crystal state, acquires a pro- 
nounced peak at zero quasi-momentum /c = 0, as shown 
in Fig. [T]3. Yet the quasi-condensation order parameter 
xf'\ namely the natural orbital of the one-body density 
matrix (OBDM) {h\^hj^cj) corresponding to the largest 
eigenvalue and hosting the condensed particles, is spa- 
tially modulated (c/. Fig. [s]^), revealing the persistence 
of solid order in the quasi-condensate. In addition, so- 
lidity can be confirmed by direct inspection of the real- 
space density (n^cr) (c/. inset of Fig. ^p). Going from 
the boundaries towards the center, the density profiles of 
both species are modulated by the crystal structure, and 
the modulation amplitudes saturate at constants which 
turn out to be independent of the system size. 

To gain further insight into the mechanism underly- 
ing the stabilization of a commensurate two-species su- 
persolid via out-of-equilibrium time evolution, we finally 
compare the equilibrium phase diagram with the non- 
equilibrium one. Fig. [2]3 shows that the superfluid/solid 
and sup erfiuid/phase- separation boundaries at equilib- 
rium overlap with the threshold of formation of the su- 
persolid out of equilibrium upon increasing This 
means that a quantum quench of the system Hamilto- 
nian to the parameter range corresponding to a super- 
fluid equilibrium ground state is a necessary condition 
for supersolidity to dynamically set in. 

The key to the dynamical emergence of a quasi- 
condensate fraction in the supersolid phase is that the 
initial trimer-crystal state (|3| has a significant overlap 
with the superfluid ground state of the final Hamilto- 
nian after the quantum quench. As shown in Fig. [4] 
for a small cluster with L = 10 sites, the ground- state 
overlap |cop remains sizable over an extended parame- 
ter range. This is intimately connected with the strong 
density-density correlations present in the equilibrium 
superfluid phase, as shown e.g. by the appearance of a 
region with super- Tonks behavior. The excellent agree- 
ment between the region featuring supersolidity and the 
region with most pronounced overlap |cop suggests the 
following mechanism: the Hamiltonian evolution follow- 
ing the quantum quench dynamically selects the ground- 
state component as the one giving the dominant contri- 



4 




FIG. 4: Overlap of the equilibrium ground state with 
the initial trimer-crystal state. The overlap |co|^ (contour 
plot) agrees well with the boundaries of the non-equilibrium 
supersolid phase (black symbols, cf. Fig. [Sj)). This suggests 
a superfiuid ground state as a necessary condition for super- 
solidity to dynamically set in. The overlap |co|^ has been 
calculated via exact diagonalization on a L = 10 chain con- 
taining three trimers. 

bution to (quasi-) long-range coherence. In essence, while 
the quantum melting phase transition occurring at equi- 
librium leads to a dichotomy between solid and superfiuid 
order, the out-of-equilibrium preparation can coherently 
admix the excited crystalline state (s) with the superfiuid 
ground state without disrupting their respective forms of 
order (see EPAPS). It is tempting to think that a similar 
preparation scheme of supersolid states can work in other 
systems displaying solid-superfluid phase boundaries at 
equilibrium. 

The supersolid transient state is an exquisitely non- 
equilibrium state, because no order can survive at finite 
temperature in ID systems with short-range interactions. 
An intriguing question arises then: in the long-time limit 
T 3fi/J| (which is only accessible numerically on very 
small system sizes) does supersolidity survive or is long- 
range order ultimately destroyed by thermalization? Re- 
cent numerical studies point towards the failure of other 
strongly correlated one-dimensional quantum systems to 
thermalize [17 . We have considered the asymptotic time 
limit using exact diagonalization for a small system (see 
EPAPS). These exact results suggest that supersolidity 



persists and the system does not converge to an equilib- 
rium thermal state (in fact even thermalization in the 
microcanonical ensemble, proposed in Ref. [18], does not 
seem to occur in our system). Whether the absence of 
thermalization survives when taking the thermodynamic 
limit remains an open question, whose answer at the mo- 
ment can only rely on experiments. 

The observation of the supersolid state prepared via 
the dynamical scheme proposed in this Letter is directly 
accessible to several setups in current optical-lattice ex- 
periments. The fundamental requirement to explore the 
phase diagrams of our model. Fig. [2j is the existence of 
a stable bosonic mixture with mass imbalance and inter- 
species interactions that can be tuned to the attractive 
regime via a Feshbach resonance. This requirement is 
met in spin mixtures of, e.g., ^^Rb atoms in different hy- 
perfine states, which acquire a spin-dependent effective 
mass when loaded in an optical lattice jT9] , and for which 
Feshbach resonances have been extensively investigated 
[20] , Moreover recently discovered Feshbach resonances 
in ultracold heteronuclear bosonic mixtures (^^Rb-^^^Cs, 
^Li-^'^Rb, ^^K-^'^Rb, ^^K-^^Rb etc. [21 , the latter re- 
cently loaded in optical lattices [22]) enlarge even fur- 
ther the number of candidate systems to implement the 
Hamiltonian ([T]). The hardcore-repulsive regime can be 
easily accessed in deep optical lattices [9] . After prepara- 
tion of the trimer crystal via an optical superlattice [15], 
the onset of coherence in the supersolid state, attained 
after a short hold time corresponding to ^ 2-3 hopping 
events of the slower particles (~ 1-10 ms), can be mon- 
itored by time-of-fiight measurements of the momentum 
distribution. The rapid onset of coherence allows the 
experimental detection of supersolidity well before deco- 
herence effects become important. On the other hand, 
the persistence of the crystalline structure can be probed 
by resonant Bragg scattering [23 . While experimentally 
the initial state will be always a mixed one and not the 
pure state in ([3| , we observe that mixedness of the initial 
state does not disrupt supersolidity in the evolved state. 

We thank J. J. Garcia-Ripoll, M. Roncaglia, and R. 
Schmied for helpful discussions. This work is supported 
by the European Union through the SCALA integrated 
project. 
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APPENDIX A: SUPPLEMENTARY MATERIAL 

In the following we present exact calculations for a 
small system which elucidate the special nature of the 
initial trimer crystal state after the quench, superimpos- 
ing the superfluid (and quasi-condensed) ground state 
with selected crystalline excited states. Furthermore, 
we compare the results for the asymptotic state of the 
time evolution with thermal states in both the canonical 
and microcanonical ensembles. Our results indicate that 
thermalization may not occur in our system. 



1. Time evolution of the initial trimer-crystal state 

We discuss here in more detail the time evolution of the 
initial trimer-crystal state into a supersolid state. The 
initial state, equation (3) of the main Letter, can be de- 
composed into the eigenstates of the final Hamiltonian 

n\Ea)=Ea\Ea) as: 



\^(t = ^))=Y,Ca\Ea). 

a 

The time-evolved state is then: 
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FIG. 5: Diagonal vs. thermal probability distribu- 
tions. The occupations of the diagonal {\ca\^ in blue) and 
canonical {\da\^ in green) ensembles are plotted as a function 
of the eigenstate energies (offset from Eqs)- Contrary to the 
thermal, continuous distribution, the trimer-crystal state em- 
phasizes certain eigenstates, while it suppresses others. The 
(superfluid) ground state contribution present in the trimer- 
crystal state is enhanced by a factor of ^ 20 compared with 
the thermal contribution. Most of the amplified excited states 
indeed show a crystalline structure with the correct periodic- 
ity, or contain density peaks at the right positions to build up 
the final crystal. Inset: The same distributions on a log-lin 
scale. The deviation of the diagonal from the thermal ensem- 
ble is even better visualized here. 



where uOa = Ea/h. 

The expectation value of any operator A can be thus 
written as 

{A), = Y.\c,\^{EM\Ea) (A3) 



^2 Re [{Ea\A\E,)clc,e'^^--^'^' 

a/6 



approaching the "diagonal ensemble" [18] or steady state 
for a large time t ^ oo. 



{A)^ = Y.\^a?{Ea\A\Ea). 



(A4) 



We now specify the discussion to the case in which the 
system is evolved with a quantum Hamiltonian whose 
ground state is both a superfluid and a quasi-condensate. 
If the initial trimer-crystal state has a significant overlap 
with the quasi-condensed ground state, namely if Cq is not 
negligible, then one can expect that the phase correlator 
of the steady state, corresponding to A = bl^bj^^^ will 
be dominated by the ground-state contribution, so that 
(quasi-) long-range order sets in. At the same time, the 
initial state has by construction a significant projection 
on excited states |£^a>o) with long-range crystalline cor- 
relations, provided that these states exist in the Hamil- 
tonian spectrum. Under this assumption, the density- 
density correlator, corresponding to ^4 = ni^cr^j,cr, will 
remain long-ranged in the steady state; this fact, com- 
bined with (quasi-) long-range phase coherence, gives rise 
to supersolidity. 

Making use of exact diagonalization on a L = 10 
chain with open boundary conditions, we have system- 
atically investigated the overlap cq between the perfect 
trimer-crystal state and the Hamiltonian ground state 



for different points in parameter space. The results are 
shown in Fig. 4 of the main Letter, and compared with 
the phase boundaries of the non-equilibrium phase di- 
agram. Fig. 2b of the main Letter. We observe that 
the non-equilibrium supersolid phase is in striking corre- 
spondence with the parameter region where cq is largest, 
suggesting that the above analysis of the onset of su- 
persolidity is quantitatively correct. Note that the time 
evolution discussed in the main Letter is restricted to fi- 
nite times, while we focus here on the asymptotic case 
t ^ oo. 



2. Comparison of the asymptotic state with 
thermal states 



The diagonal-ensemble expectation value of equation 



(A4) is here compared with a thermal average in the 
canonical ensemble 



{Ah 



E 



\da\^{Ea\A\Ea) 



(A5) 



with I dap = ex.p{—Ea/kBT)/Z the Boltzmann weights, 
ks Boltzmann's constant, T the temperature and Z = 
'^^ex.-p{—Ea/kBT) the normalizing partition function. 
In addition, we introduce for comparison the statistical 
average in the microcanonical ensemble 

{A)E,„,dE= J2 l/N^{EM\Ea), (A6) 

which averages over eigenstates within an energy window 
±dE around the inital energy E^^ = {"^{t = 0)\H\'^{t = 
0)). Njn is the number of eigenstates contained in that 
energy window. 
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In order to compare the diagonal with the canonical 
and microcanonical ensembles, we have chosen to exactly 
diagonalize a system of three trimers (N^ = 6^N^ = 3) 
in an open chain of L = 10 sites. We present in the fol- 
lowing results for the parameter pair (J| = 0.2J|,[/ = 
3J|), where super solidity exists according to our non- 
equilibrium phase diagram, Fig. 2b of the main Letter. 
Under these conditions, the ground state energy yields 
Egs — — 13.2J|, while the initial trimer-crystal state car- 
ries an energy E^^ = — 12J|. In order to determine the 
correct temperature for the canonical ensemble, we have 
varied T until the condition (H)t = ^in was met. This 
analysis yielded ksT :^ 0.82J|, which we use henceforth 
for the comparison with the canonical averages. 
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FIG. 6: Real-space density (nj) in all three ensembles. 

While the diagonal ensemble (n j ) oo (blue) shows a clear crys- 
talline pattern, this structure is washed out completely in the 
canonical ensemble {nl)T^o.82j^ /ks (green). Results for the 
microcanonical ensemble {n\)Ei^,dE are shown for energy win- 
dows dE — 0.2J| (red) and dE — 0.6J| (cyan). All thermal 
ensembles deviate strongly from the density structure at time 
t ^ oo (diagonal ensemble). 
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FIG. 7: Momentum profile {ni) in all three ensembles. 

Due to the significant weight attributed to the ground state, 
the diagonal ensemble {n\^)oo (blue) features an enhanced 
quasi-condensation peak at /c = 0. This peak is suppressed in 
all thermal ensembles {n\^)T=Q.s2j^ /ks {ni)Eo,dE (same 
colouring scheme as in Fig. [g]). 

Fig. [5] compares the diagonal ensemble induced by 
the initial trimer-crystal state with a thermal, canonical 
ensemble. The trimer-crystal state has a finite projec- 
tion on the quasi-condensed ground state as well as on 
distinct excited states. Further inspection into those ex- 
cited states shows that their characteristic density pro- 
files matches the crystal structure of the initial state. 
Hence the selection of excited states in the diagonal en- 
semble is fundamentally governed by the broken transla- 
tional invariance present in the initial state. In contrast. 



the canonical ensemble averages over all eigenstates re- 
gardless of their displaying crystalline order, a fact which 
makes the loss of the crystalline structure unavoidable. 

The density profiles (for the |-bosons) shown in Fig. 
|6] corroborate the previous statements. The diagonal en- 
semble induced by the trimer crystal is compared here 
with thermal averages in both the canonical and micro- 
canonical ensembles. While the density profile in the di- 
agonal ensemble still displays the "memory effect" of the 
initial crystalline state, the thermal states exhibit only 
small density modulations (in the microcanonical ensem- 
ble) or no modulation at all (in the canonical ensemble). 
Furthermore, the momentum profiles shown in Fig. [7| un- 
derline a non-thermalization scheme of the time-evolved 
crystal state. While the density profile of the diagonal 
ensemble exhibits a pronounced peak at quasimomentum 
/c = 0, this peak is almost completely washed out for the 
thermal ensembles. 

In view of the two observables discussed here, a 
thermalization of the evolved trimer crystal state can 
be excluded, at least for the finite-size system we are 
considering. This confirms the observations of "non- 
thermalization" in other ID finite-size systems [17]. 

Our exact diagonalization study is limited to a small 
cluster, and it cannot exclude a priori that thermalization 
appears for larger system sizes: this would require that 
the diagonal ensemble converges to the microcanonical 
one, which ultimately converges to the canonical ensem- 
ble in the thermodynamic limit. 

3. Numerical results of long-time evolutions 

Here we present an example of the scaling analysis for 
the results of a long-time evolution up to r = 150fi/J|. 
Fig. |8] shows that observables averaged over the last 
r/4 interval of the time evolution display the charac- 
teristic one-dimensional supersolid scaling, analogous to 
- but much more marked than - the one observed at 
short times (compare Fig. 3a-b of the main Letter). In- 
deed we observe a linear scaling of the structure factor 
peak S{qtr = 27r/3) with system size, typical of solid or- 
der, and an algebraic sub-linear scaling of the condensed 
atoms, signaling quasi-condensation. Repeating this scal- 
ing analysis for a fine mesh of parameter space leads to 
the confirmation of the supersolid phase shown in Fig. 
2b of the main Letter. 

A word of caution is necessary in the case of long-time 
evolutions. The truncation of the Hilbert space, inher- 
ent in all numerical algorithms for time evolution not 
based on full exact diagonalization [16], has the general 
effect that the accuracy of the results rapidly degrades 
with time, and the instantaneous measurements become 
practically unreliable in the long-time limit. We observe, 
however, that observables averaged over time intervals 
> r/10 do converge with high precision upon variation 
of the bond dimension D. These time averaged results 
are indeed the object of the scaling analysis in Fig. [5] 
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FIG. 8: Scaling analysis of the long-time evolution 

data, (a) Structure factor peak S{qtr = 27r/3); (b) Quasi- 
condensate peak (nfc=o)- Boxes (diamonds) stand for particle 
species | (t), respectively. Parameters: Jj/Jt — 0.15, U/J^ = 
2.5 (blue symbols) and JJJ^ = 0.40, C// Jt = 9.0 (red sym- 
bols). 



